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Abstract: We describe an algorithm for improving subsequent parton shower emis- 
sions by full SU(3) color correlations in the framework of a dipole-type shower. As 
a proof of concept, we present results from the first implementation of such an algo- 
rithm for a final state shower. The corrections are found to be small for event shapes 



and jet rates but can be more significant for tailored observables. 
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1. Introduction 

Parton showers and event generators are indispensable tools for predicting and un- 
derstanding collider results [1-3]. Considering their importance for interpreting LHC 
results, it is essential to have a good understanding of their approximations and lim- 
itations. 

Significant progress has been made over the last years in the areas of matrix 
element merging at leading order [4-9], and matching at next to leading order [10-15], 
as well as improvements to shower algorithms using dipole- type evolution [16-20], 
partly generalizing ideas presented originally in [21]. There has also been theoretical 
progress in the direction of treating full SU(3) at the level of shower evolution [22], 
and at the level of dealing with the color space [22-24]. 

In dipole-type parton showers, only large N c color connected partons 1 can radiate 
coherently. The, typically l/Af? suppressed, coherent emission from other non-color 
neighboring partons is neglected. However, as the number of perturbatively emitted 
partons n increases, the number of possible contributions from non color neighbors 
grows roughly like n 2 /2, while the number of possible coherent emissions from color 
neighbors grows only like n. With the higher multiplicities at high energy hadron 
colliders, these color suppressed contributions may thus become important, though 
this counting does not include the dynamics of multiple parton emission. A treatment 

11 Color connected' here means that two partons share a common color line as used in the double 
line notation. Note that the double line notation can be extended beyond the large N c limit when 
including an appropriate singlet contribution for the color structure of the each gluon propagator. 
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of these color suppressed contributions within a full-fledged shower framework is thus 
necessary. 

This paper reports on a first step in the direction of incorporating color sup- 
pressed contributions by iterative use of color matrix element corrections. We outline 
an algorithm to calculate these color matrix element corrections to improve shower 
emissions by the full SU(3) color correlations and present a proof of concept imple- 
mentation based on the shower implementation described in [15,20]. The first results 
are presented for e + e _ — > jets. 

For the implementation, a C++ color algebra package, ColorFull [25], has been 
developed for keeping track of the color structure and aiding calculation of high 
multiplicity matrix elements squared and color correlated matrix elements. It has 
also been necessary to consider the treatment of negative weights in the context of 
partons showers, as there are negative contributions to the radiation probability from 
color suppressed terms [26]. 

This paper is organized as follows: in Sec. ^| we recapitulate the dipole factor- 
ization scheme and the dipole shower evolution. In Sec. [3] we define an algorithm 
for the inclusion of full color correlations for subsequent parton shower emissions. 
The details of the color space treatment are given in Sec. |j. In Sec. we present 
numerical results from the implementation of a parton shower improved by full color 
correlations. Sec. || draws some conclusions and gives an outlook on future develop- 
ments. 



2. Dipole Factorization and Dipole Shower Evolution 

Dipole factorization, [27,28], states that the behavior of QCD tree- level matrix ele- 
ments squared in any singly unresolved limit involving two partons i,j (i.e. whenever 
i and j become collinear or one of them soft), can be cast into the form 

\M n +i(-~,Pi, -,Pj, ...,Pk, ■■■)\ 2 ~ 

Y] 7T~ ( M n{---,Pij, ■■■iPh ■■■)\ V ijAPi,PjiPk)\Mn{...,p i ~ j , ...,p~ k , ...)) , (2.1) 

tr" . -i>i ■ Pj 

where \M n ) - which is a vector in the space of helicity and color configurations - 
denotes the amplitude for an n-parton final state. Here an emitter ij undergoes 
splitting to two partons i and j in the presence of a spectator k which absorbs the 
longitudinal recoil of the splitting, k —> k, such as to keep the momenta both before 
and after emission on their mass shell while maintaining exact energy-momentum 
conservation. Considering gluon emission only, this can be interpreted as a dipole 
ij,k emitting a gluon of momentum pj. (However, g — )■ qq splittings fit into this 
framework without conceptional changes). 
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Figure 1: Sample diagrams contributing to the dipole kernels Vy jt and Vjyi in the 
singular limits for gluon emission off a quark or antiquark. Similar diagrams are present 
for gluon splittings. 

In general, the insertion operator V contains both color correlations stemming 
from soft gluon emissions and spin correlations originating in the collinear splitting 
of a gluon. We shall limit ourselves to massless quarks and color correlations only, 
using the spin averaged dipole splitting functions presented in [28]. In this case 



T 2 -. 



(2.2) 



in terms of the color charge operators Tj 2 . This notation is independent of the basis 
considered for color space, though we shall stick to one particular choice of basis, to 
be discussed in Sec. f|. 

The kernel Vy^ only contains a collinear enhancement with respect to the split- 
ting ij — > i,j. It is constructed by rearranging singular contributions of sets of 
diagrams as depicted in Fig. [I] while making use of the fact that the amplitude is a 
color singlet, 



rj-i2 



(2.3) 



For gluon emission off a final-final dipole configuration, the case considered in the 
present paper, the kernels read 



Vqg,k(Pi,Pj,Pk) = C f 

Vgg tk (Pi,Pj,Pk) = 2C A 



2(1 - Z) 



z) 2 +p]_/s ijk 
1-z 



[1 + z) 



1 



+ p 2 ± /s 



+ 



ijk 



z 2 +p]_/s ijk 



2 + z(l 



(2.4) 
:) 1(2.5) 



where we have introduced the relative scalar transverse momentum p± and longitu- 
dinal momentum fraction z of parton i as given by the Sudakov decomposition 



2 We use the conventions on the color charge algebra as given in [28]. Note that the Casimir 
operators T 2 = Cf if i is an (anti-)quark, and T 2 = Ca if « is a gluon are included in the definition 
of the splitting kernels V; thus the color correlations are normalized accordingly. This convention 
is more transparent when comparing to parton showers in the large N c limit to be discussed below. 
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Pi = zpq + —Pi + k± (2.6) 
Pi = (1 " *)PQ + ,*\ Pk ~ kx (2.7) 

"=( X -5(T^)«- (2 ' 8) 

with p 2 r. = p~ = 0, a space like transverse momentum k± with k\ = —p\ and k\-p^ = 
k±-pj t = 0. With this parametrization we also have = (Pi+Pj+Pk) 2 = (P{j+Pk) 2 - 
Reinterpreting the dipole factorization in terms of emitting dipoles, parton cas- 
cades can be based on it and, after slight modifications to initial state radiation, have 
been shown to exhibit the proper features in terms of coherent gluon emission and 
logarithmic accuracy [20]. The present work is an extension to the implementation 
of a coherent dipole evolution, details of which have been reported in [15]. Besides 
using spin averaged dipole kernels - like any other parton shower implementation 
available so far - the Monte Carlo implementation in [15] treats the color correlation 
operator for the shower in the large N c limit i.e. when using the double line notation 
of color 



TV. • T fc 1 

— > —5(ij,k color connected) , (2.9) 



T 2 , ' l + 5r 

where 8^ = 1 if ij is a gluon, and vanishes otherwise. Here the initial assignment and 
further evolution of large N c color charge flow is extensively discussed in e.g. [15]. 
Let us only recall here that to each of the two emitter-spectator partitions of a large 
iV c color connected pair of partons we connect an emission probability related to the 
dipole kernels as 



dPij,k(p±, z;Pij,P~ k ) = d*J(p±, z;Pij,Pk) v ij,k(pl, z 'iPivP~k) ■ ( 2 - 10 ) 

Z7T p . 



Here, J represents the Jacobian for changing variables from Pi,Pj,Pk to p-j, p^, p±, 
z and an azimuthal orientation, which has been integrated over. For a final-final 
dipole configuration, J(p±,z;pr.,p~) = (1 -p]_/z(l - z)s ijk ). 

In this paper we report on generalizing this picture from large N c color connected 
dipoles to all pairs of partons, taking into account the exact evaluation of the color 
correlations. Thus, the notion of dipole chains and their splittings does not apply 
anymore. Instead every pair of partons has to be considered competing for the next 
emission. Details of the generalized algorithm are given in Sec. H below. 
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3. Color Matrix Element Corrections 



Following the arguments of how a parton cascade is constructed from the dipole fac- 
torization formula fl2.1|) , we may recast fl2.10| ) to include the exact color correlations 
by assigning splitting rates 

%f*J<A,r, n) ,riV w .< J l,r, n3 ,ri x ^^0^ (3.1) 
to each pair ij, k of partons to generate an emission j off an n-parton system. In 



the large N c limit the factor after the times sign reduces to (|2.9|) for color connected 
pairs of partons and zero otherwise. Here, by allowing any pair of partons to radiate, 
it accounts for correcting the emission to include the exact color correlations. Owing 
to the similarity to so-called matrix element corrections present in parton showers, 
[29,30], we refer to this improvement as color matrix element corrections. 

At the level of the hard process, from which the shower evolution is starting, 
the calculation of \Ai n ) is straightforward, but once the first emission is performed, 
the issue of how to describe the evolved state has to be addressed. By picking a 
particular basis for the color structures, we can map \JA n ) to a complex vector 
Ai n G C dn , where d n is the dimensionality of the basis for n partons, 

\M n ) = ^c n , a |a n ) M n = (c n ,i, ...,c n4n ) T . (3.2) 

a=l 

Then we have 

\M n \ 2 = MiS n M n = Tr (S n x M n Mi) (3.3) 

with S n being the scalar product matrix, (S n ) a p = (a n \f3 n ), for color basis vectors 
\a n ) and \/3 n ). Moreover, the color correlated matrix element for emission from ij 
and k can be written as 

(M n \T {j ■ T- k \M n ) = Tr (S n+1 x T kn M n MlTt n ) (3.4) 

in terms of matrix representations of Tj--,Tj. 6 C dn+1 ' dn . These matrices describe 
the color space effect of emission from the partons ij and k by mapping a basis 
tensor in n-parton space to a linear combination of basis tensors in n + 1-parton 
space. In the "trace basis", which is used here, (c.f. Sec. |]) the TV- matrices are 
very sparse [22,23]. The scalar product matrices, S n , S n+ i, are however dense, and 
the calculations of these matrices is a key ingredient for running a parton shower 
improved by color matrix element corrections. 

As we want to keep the full color structure the evolution is (analogous to [22]) 
done keeping the amplitude information, using a matrix M n , which initially - for 
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the hard process the shower starts with - is given by M n = M. n Jv[\. After having 
performed an emission with momentum pj off an n-parton configuration, this quan- 
tity for the n + 1 parton configuration can then be obtained from the spin-averaged 
dipole kernel as 

^ = -££ f^liMl T- kn M n Tj . (3.5) 

Matrix elements squared and color correlated matrix elements for each subsequent 
color matrix element correction are then calculated as in ( |3.3|) and (|3.4j ) upon re- 
placing AinM^ — > M n . One property of the color matrix element correction weight 
is that it is not necessarily positive definite for subleading-iV c contributions. The 
authors have shown in [26], that this does not pose a problem for a Monte Carlo 
implementation. In particular, from the very definition of the dipole factorization 
as an approximation to a tree-level matrix element squared, the sum of the splitting 
probabilities including all pairs i, k is strictly positive definite. Thus we can imme- 
diately apply the interleaved competition/ veto algorithm [26] to generate kinematic 
variables with the desired density. 

More precisely we select a set of candidate emissions ij,k — > i,j,k at scales 
P±tj,k ( wnen starting from a scale Q\ associated with the dipole ij, k) according to 
the Sudakov form factor 

-\nA ij!k (p Xijk \Q x ) = / — / dzV ijtk (ql,z;pr p~ k ), (3.6) 

2vr J v 2 q, J z _t^,nn.\ 



. h q± Jz_(ci,Qi) 
if Vij s k is positive. Here, in accordance with ( |3.1| ), 

(3-7) 

and z±(p 2 L ,Q 2 L ) = (1 ± a/1 — p 2 L /Q 2 L )/2. Amongst the candidate splittings of ij,k, 
we pick the one with largest p 2 .-. ,= p\ and accept the corresponding splitting with 
probability 

J2{j J2k^ij ^ij,k {p\ j z \ Pij :Pk) ^ 



Eij Ek^j P*i,fc(pL z\ Pi^PkWi^jAph z\ Pij, Pi)) ' 

Upon rejection, the selection is repeated, setting Q± = p±, until the next p± is 
accepted or eventually p± is found at the infrared cutoff \i = (9(1 GeV). Upon 
accepting the candidate splitting with scale p±, the dipole ij,k is chosen to define 
the Sudakov decomposition of the emission kinematics. The generated emission is 
then inserted into the event record, and M n+ \ is determined using the generated 
kinematics according to (O). 
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Figure 2: Examples of color structures, (a): An open quark line with 6 gluons cor- 
responding to the color structure t^tq^tqlq^q^tqlq^ 9 ^. This is the only type of color 
structure needed for e + e~ — >• qq — >• qq + N g gluons at leading order, (b) A closed quark 
line corresponding to the color structure t q ^ qi tf 1 q 2 tf 2 q 3 t q 'tqi4lq5 t qtq6 = Tr [t 9l t 92 t 9 H 9 H 9 n 9fi }. 

4. Color Basis Treatment 

Throughout this paper we use the basis obtained by first connecting all gg-pairs in 
all possible ways (for one gg-pair only one way), and then attaching gluons to the 
quark-lines in all possible ways (N g \ ways for one gg-pair), see Fig. |a. This basis is 
sufficient as long as only leading order QCD processes are considered. In the presence 
of (QCD) virtual corrections, or processes mediated by electroweak (or other color 
singlet) exchanges, we would also have to consider basis vectors obtained from direct 
products of open and closed quark-lines, as depicted in Fig. |[ 
Due to the relation 

= ^(Tr[rt c ] - Tr[tW]) = " ^ 

J-R J-R 

where Tr = Tr[t a t a ] (no sum) is taken to be 1/2, the effect of gluon emission and 
gluon exchange is trivial in this basis [22,23]. Graphically the emission process can 
be exemplified as in Fig. |[ Using the notation 

\{qgi92-g m q}) = CC-Ca ( 42 ) 
we may in general write for the emission of gluon g n +\ from the gluon at place % 

\{qg 1 ...g i ...g m q}) -»■ \{q g 1 ...g i g m+1 ...g m q}} - \{qgi...g m +i 9i-9m q}) (4.3) 

where the overall sign depends on the sign convention for the triple gluon vertex, 
and thus has to be matched with the convention for the momentum dependent part. 
Here we choose the triple gluon vertex factor to be ifg, g ,g m+1 where g denotes the 
emitters color index before emission, g denotes the emitters color after emission, and 
g m+ i denotes the color index of the emitted gluon. For quarks and anti-quarks the 
effects of gluon emission are similarly 

\{qgi-g m q}) ->■ |{?3m+i3i-fa«}> forg (4.4) 
\{qgi...g m q}) -> -\{qgi-9m9m+iq}) forg. (4.5) 



g 9 g q g 1 g 9 g 4 g Q The disadvantages of this type 

of basis - obtained by connecting 
partons in all possible ways - is that 





9 1 9 2 9 4 93 9i 94 9 2 93 [t is orth °g° nal onl y in the N c oo 

limit, and that it is overcomplete 
for finite N c . For N c = 3 the low- 
est number of partons for which this 
Figure 3: The graphical representation of ( p~3| ) kind of basis is overcomplete is four 
for i = 2. gluons (incoming plus outgoing). 

We note that there are two com- 
putational steps which naively scale as (N g \) 2 . The first step is the step of calculating 
all the scalar products between the N g \ basis vectors. The second is the step of iden- 
tifying the new basis vectors when a given parton in a given basis vector has emitted 
a new gluon; the number of basis vectors in the initial basis scales as (N g — 1)! (start- 
ing from N g — 1 gluons), and the number of basis vectors in the final basis as N g \. 
This would give an overall scaling behavior of (N g — l)\(N g + l)N g \ ~ (N g \) 2 , as there 
are (N g + 1) possible emitters. 

However, enumerating the basis states in a unique way, it is possible to calculate 
the numbers of the new basis vectors when a gluon has been emitted from a given 
parton in a given basis vector. As there is no need to compare the result after 
emission to all basis vectors in the iV 9 -basis, this reduces the computational effort to 
scale rather as (N g — l)!(iV 9 + 1) in the step of identifying all new color states when 
starting in (N g — 1)! possible basis vectors and emitting from N g + 1 possible emitters. 
This is clearly much better than the {N g \) 2 scaling for calculating the scalar products, 
which thus is the major bottleneck for sufficiently many partons. Fortunately these 
calculations need only be performed once, and can then be stored numerically for 
use in the color matrix element corrections. 

The scalar product matrices have been checked to agree with scalar product 
matrices calculated by the Mathematica code published along with [24] for up to 
4 gluons, and against a new (yet unpublished) Mathematica package. Similarly the 
matrices representing TV- • T^. have been checked to agree for up to 3 gluons compared 
to the old Mathematica code, and for up to 4 gluons compared to the new. 



5. Results 

We here discuss first results from a subleading N c improved final state shower, origi- 
nating from e + e~ — > qq at 91 GeV center of mass energy. We consider gluon emission 
only, as there is no soft enhancement present for a g — > qq splitting. The strong cou- 
pling constant is taken to be fixed, a s = 0.112. The subsequent gluon emissions are 
performed as described in Sec. [3] using one of three options for the color structure: 
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1. Full color structure; the full SU(3) color structure with splitting kernels as in 
(|3.1|) is used with emissions generated as in Sec. [| 



2. Shower color; the color matrix element correction is evaluated in the large N c 
limit, though the color factor Cp entering q —)■ qg splittings is kept at its exact 
value Cp = 4/3. This resembles current parton shower implementations. 

3. Strict large iV c , all iV c suppressed terms are dropped, implying Cp = 3/2. 

Before turning to the numerical results it is instructive to study the expected effects 
already at an analytical level, and to clarify how the shower limit maps to the stan- 
dard shower implementation. We therefore consider the correction weights for up to 
three emissions, labelling each intermediate state with N g = n — 2 additional gluons 
as qiq29z---9n- The gluon g n , is always considered to be the one emitted in the last 
transition, and the weight associated to a dipole i, k in an ensemble with n partons 
is denoted by 

s V in>k (pi,p n ,p k ) = V£ . (5.1) 



The color matrix element correction for emission off a dipole i,k in an n-parton 



ensemble is denoted 

1 (M n \Ti ■ T k \M n ) 1 Tr (^n+i x T ^ M " T *U) 
"T? \Mj 2 = ~Tf Tr(S n xM n ) 



< (5-2) 



for brevity. Keeping the full correlations we find the corrections for the first two 
emissions to be given by 



«>12 = W 21 = 1 ( 5 -3) 

n n 9 



w 13 = w 23 



w 31 = w 32 



Wf 2 = W\y 



The negative contribution from the qq dipole in the qqg system has already been 
noted in [31]. Note that, e.g. wf 2 + wf 3 = 1 as dictated by color conservation. In 
the shower approximation we have 

w 2 i2 = w 2 2l = l (5.4) 

^13 = wl 3 = l 



w 31 = w 32 



W 12 = ^21 = 
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matching precisely the naive expectations on the subleading iV c contributions, that 
there is no radiation off the qq dipole in a qqg system. For the first emission there is no 
difference between the two approximations reflecting the triviality of the color basis 
in that case. Also, gluon splittings in a qqg system do not exhibit any subleading N c 
correction. 

More non-trivial dynamics are present for the fourth emission. Introducing the 
relative magnitudes of dipole kernels encountered in the four parton system with 
respect to the q(q)g dipoles, 

rik = v? 3 + v 3 \ + vg + vl 

we find for the exact correlations 

W t 3 = W 24 = 7 1 — ( 5 - 6 ) 



W 24 = 
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r 23 - 


-\ri2 
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1 - 


\r l2 


W 23 = 
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ri3 ~ 


-\r 12 
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1 - 


\r l2 


w\ 2 = 
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r 23 - 


-\r 12 
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1 - 


\r l2 


™32 = 


1 


ri3 ~ 


-\r l2 


2 


1 - 


\r l2 



w u = 



U>31 



w 41 = 



""34 = -4, = J 

4 4 1 1 - f r i2 

Win = U> 91 = T 

12 21 8 1 - \r 12 



We note that the correction weights only depend on the quantities r^. Conversely, 
in the large N c limit we find 

(5.7) 



«4 


= «4 


= r 23 




= W 23 


= ri 3 




= «4 


1 

= 2 r23 






1 

= 2 ri3 


«& 


= «4 


1 

~~ 2 




= W 21 


= . 



Again, there is no radiation off the qq dipole, the gg dipole precisely matches the 
standard shower implementation. The different combinations of q(q)g dipoles resem- 
ble the standard shower implementation only in the case that - in the three parton 
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system - the splitting functions of one dipole had been much larger than the split- 
ting functions of the other dipole. In this case the weights precisely match up the 
distribution of radiation generated by the shower once it has decided which dipole 
was to radiate - this history is of course closely linked to the hierarchy of splitting 
kernel values encountered. Indeed, the deviations between the shower approximation 
and the standard shower implementation have been found to be negligible. 

The shower is terminated when either the iV ff :th gluon (we here consider up to 
N g = 4, 5, 6 gluon emissions) is emitted, or when the infrared cut-off, taken to be 
1 GeV, is reached. Only about 1% of all events radiate up to 6 gluons. 

As, due to computational limita- 
tions, we only consider a limited num- 
ber of emissions, we have, for the 
observables considered here, checked 
that the prediction stabilizes when go- 
ing up to six emissions. The distri- 
bution of the number of emissions is 
shown in Fig. |] for the various approx- 
imations considered. We find that 
taking the strict large N c limit rises 
the probability for many gluon emis- 
sions by more than 40%. This can 
in part be attributed to the fact that 
C F = T R (N 2 C - l)/N c , which equals 
4/3 (using T R = 1/2) is replaced by 
its leading part 3/2 also for collinear 
emissions. The buildup of this in- 
creased radiation probability accounts entirely for the difference between the "strict 
large N c " and the "shower" treatments. Apart from this effect, there is another effect 
which makes both the strict large N c and the shower treatments of color, result in 
more radiation, compared to keeping the full color structure. This may be attributed 
to the fact that the interference between color structures, where two gluons cross each 
other, as depicted in in Fig. |5](b), comes with negative sign. These terms - which 
are not present in standard parton showers - thus seem to lower the probability for 
radiation. For a large number of emissions, color structures where gluons cross 
each other in more complicated ways dominate the color suppressed contributions. 
Starting from two gg-pairs there are terms which are only suppressed by one power of 
N c . We thus caution that the effects of color suppressed terms may be significantly 
larger at the LHC. 

For the LEP-like setting considered here, we have investigated a set of stan- 
dard observables, event shapes and jet rates. In all cases we find the deviations of 
the shower approximation to be small, up to a few percent, when compared to the 
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Figure 4: The normalized distribution of the 
number of emissions above the infrared cutoff 
comparing the different levels of approximation. 
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(a) (b) (c) 

Figure 5: With diagrammatic techniques an amplitude square, as in (a) can easily be 
evaluated to N c Cp 9 . The type of interference depicted in (b), where 2 gluons cross at 
one point contributes negatively to the radiation probability with —Tp>C F 9 1 , whereas the 
interference depicted in (c) contributes positively with T^Cp 9 2 (N^ + 1)/N . 
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Figure 6: The distribution of one minus thrust, r = 1 — T. The left panel compares 
different approximations of the color structure. The right panel compares predictions de- 
pending on the number of emissions using the full color structure. The prediction stabilizes 
from five emissions onward, as is the case for the other observables considered. 



full color treatment. The differences between the strict large N c and the "shower" 
treatment option are often larger, in the cases of thrust, T, and sphericity up to 
roughly 10%. In Fig. |B] we show the r = 1 — T distributions for the three options 
of color structure (left), and - using the full color structure - depending on how 
many emissions we allow (right). The right plot shows that the prediction stabilizes 
from five emissions onward. We note that there is not much difference between the 
full and the shower-like treatment of the color structure, whereas the strict large N c 
treatment results in less pencil like events. This can be understood by noting that 
there are more emissions in this indicated in Fig. |J 

We now turn our attention to differential Durham n-jet rates, Fig. [7], showing 
the distribution of the resolutions scale, y = 2 mm(Ef, Ej)(l — cos0y)/s, at which 
an n+ 1-jet event changes to an n-jet event. We note that for the two-jet rate, there 
is almost no difference between the full and "shower" option whereas, increasing 
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Figure 7: Durham differential n-jet rates for transition from three to two (left), and six 
to five jets (right). 

the number of jets, we may see a few percent difference, although this difference 
may be just a statistical fluctuation. For the two jet rate the small difference can 
be understood by considering that the transition from a three to a two-jet event 
is only sensitive to the emission of a gluon from a quark-line, i.e. to Cf, which is 
correctly described in standard showers, as well as for the case that no non-trivial 
color correlations are present for the first emission. Nevertheless, the difference 
between shower and full remains very small, at most a few percent, for the differential 
n-jet observables. 

In general we find that the differences between the shower approximation and full 
correlations are small for event shapes and jet rates considered here. The explanation 
for this is likely that these observables either are mainly sensitive to the collinear 
singularity - which is treated with the correct color factor in standard parton showers 
- or are mainly sensitive to the first hard emission, in case of which there are no non- 
trivial correlations present. 

The next candidates to check for larger effects are four-jet correlations, which 
have been studied at LEP, mainly to investigate the non-abelian nature of QCD, i.e. 
these are all very sensitive to g — > gg splittings. An example, the distribution of the 
cosine of the angle between the softest two jets in four-jet events at a Durham-jet 
resolution of y = 0.008 are shown in Fig. |||. No large deviations are observed between 
the different approximations. A closer consideration of the color space for a qq pair 
and two gluons reveals that this may actually be expected. Note that there is almost 
no difference between the shower and strict large N c approximations, which can be 
attributed to the fact that these observables mainly probe gluon splitting which is 
not sensitive to the difference between these approximations. 
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Figure 9: The average transverse momentum and rapidity of partons four onward w.r.t. 
the thrust axis defined by the three hardest partons. 
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To investigate the effects of soft 
coherent emissions, we have therefore 
also studied the average transverse mo- 
mentum and rapidity of parton four 
onward with respect to the thrust axis 
defined by the three hardest partons. 
The result is shown in Fig. ^. Here 
the effects are much larger. We find 
that the average transverse momentum 
is harder if the full color structure is 
kept, for large (p±) up to 20%. We also 
see that the strict large N c approxima- 
tion is somewhat closer to the predic- 
tion including the full correlations. We 
leave a more detailed study of this class 
of observables as future work, includ- 
ing an operational definition in terms 
of anti-/c_i_-type jets to build up a reference system from hard, collinear jets while 
being able to look at the orientation of soft radiation relative to this system. 




-0.5 





COS Q;m 
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Figure 8: Distribution of the cosine of the 
angle between the softest jets in four-jet events. 



6. Conclusions and Outlook 

In this paper the first results from a subleading N c parton shower were presented. We 
have considered final state gluon emissions using iterative matrix element corrections 
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to treat the full SU(3) structure for each subsequent emission. 

A major conclusion is that standard LEP-observables, including event shapes 
and jet rates, are only affected by at most a few percent. For a class of observables 
dividing the event into a hard reference system and accompanying radiation we 
expect larger differences. This is indeed seen for tailored observables like the average 
transverse momentum with respect to a thrust axis determined by a system of three 
hardest partons, in case of which we see deviations as large as 20%. 

The small differences can, in the case of the standard observables, largely be 
explained by the fact that these observables are either sensitive to collinear radiation, 
or to the first hard emissions. 

Another contributing factor is that the color suppressed terms can be seen to be 
quite small from considerations of color space alone, when starting from a gg-pair. 
This is not the case if the hard scattering process is e.g. QCD 2 -> 2. In this case we 
may see more striking differences between the approximations considered, though we 
cannot yet make a definite statement. The simulation framework is general enough 
to cope with this case and we leave a detailed discussion of subleading N c effects at 
hadron colliders to future work. 

In addition to what is presented here, there are several other effects which should 
be included before it can be claimed that a parton shower fully simulates SU(3) 
physics. Apart from including a running coupling constant and the g — >■ qq split- 
ting kernel, we like to include virtual color rearranging gluon exchanges. We view 
this work as a proof of concept, and a first step towards quantifying the impact of 
subleading N c contributions. 

Finally, the parton shower outcome should be hadronized. While this is in itself 
an interesting task, it is much beyond the scope of the present paper. In an ordinary 
parton shower, where the shower outcome corresponds to a well defined probabilistic 
color line arrangement, this color structure is fed into the hadronization model. 
Here, the state after showering contains amplitude level information, and can thus 
not simply be input into existing hadronization models. While studying the influence 
of the hadronization on the parton shower outcome is an interesting task, we thus 
refrain from it at this stage. 
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